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SOUND  PROPAGATION  IN  RAREFIED  OASES 


ABSTRACT 

A  study  of  sound  propagation  in  rarefied  gases  Is  described  with 
particular  emphasis  on  the  propagation  In  the  geometrical  relaxa¬ 
tion  regime  where  the  sound  carrying  molecules  do  not  suffer 
lntermolecular  collisions  during  the  flight  from  the  transmitter 
to  the  receiver.  The  propagation  Is  characterized  by  the  atten¬ 
uation  and  the  phase  parameters.  The  dependence  of  these 
dispersion  parameters  on  the  separation  between  the  transducers, 
the  velocity  distribution  of  the  molecules  subsequent  to  their 
Interaction  with  the  transmitter  surface  and  the  boundary  condi¬ 
tions  describing  molecular -surface  Interactions  Is  ascertained  In 
terms  of  a  semi -phenomenological  theory.  Experiments  for  measuring 
the  dispersion  parameters  are  described  and  the  results  are  dis¬ 
cussed  and  compared  with  theory.  A  theory  and  an  experiment  for 
determining  surface  accommodation  coefficients  using  sound  propa¬ 
gation  In  rarefied  gas  are  outlined.  Finally  a  theoretical  model 
Is  constructed  that  leads  to  a  reasonable  qualitative  and  quanti¬ 
tative  description  of  sound  propagation  In  rarefied  gases. 
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I.  INTRODUCTION 

The  propagation  of  sound  In  gases,  especially  In  rarefied  gases, 
offers  a  method  for  studying  gas  dynamics  on  a  molecular  basis. 

Prom  the  theoretical  point  of  view  this  study  Is  particularly 
attractive  since  the  amplitudes  which  define  the  sound  field  are 
usually  small  compared  with  the  amplitudes  which  define  the  basic 
motion  of  the  molecules.  Thus,  one  is  able  to  linearize  the  kinetic 
equation  of  motion.  Of  particular  Importance  Is  the  linearization 
of  the  Boltzmann  equation 


(&  +  i)  “(**i*t>  ’ 


(1.1) 


I  is  the  source  tern  arising  from  inte molecular  interactions, 
mm 

1^  is  the  source  term  arising  from  molecular-surface  interaction, 
and  i  is  the  molecular  velocity  vector.  The  function  h  represents 
the  perturbation  tern  in  the  distribution  function  which  is  caused 
by  the  superposition  of  the  sound  field. 


A  considerable  volume  of  literature-^/  exists  which  deals  with 
methods  and  techniques  for  obtaining  solutions  to  the  linearized 
Boltzmann  equation.  These  analyses  are  confined  chiefly  to  the 
' j. r  field  where  all  memory  of  surface  interaction  has  been  de¬ 
stroyed  by  intennolecular  encounters.  One  considers  a  source 
tern  in  the  Boltzmann  equation,  devoid  of  any  molecular-surface 
interaction  contribution  (1^  ■  0).  The  results  which  are  obtained 

indicate  that  measurements  of  sound  in  the  far  field  are  sensitive 
to  some  of  the  properties  of  intennolecular  collisions. 
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The  analysis  cf  sound  propagation  sufficiently  near  a  transmitter 
surface  such  that  molecules  which  carry  the  sound  from  that  sur¬ 
face  arrive  at  the  observation  point  without  undergoing  inter- 
molecular  collisions  can  also  be  performed.  In  this  analysis  the 
intermolecular  encounters  may  be  ignored  (I  -  0).  Consequently 

one  expects  the  solution  to  the  equation  of  motion  to  be  dominated 
by  molecular-surface  interaction  parameters.  An  experimental 
study  of  molecular-surface  interaction  by  acoustical  means  is  thus 
suggested. 

In  considering  the  kind  of  information  which  one  may  reasonably 
expect  to  obtain  from  experiments,  it  is  convenient  to  consider 
the  properties  of  sound  propagation  in  three  distinct  regimes  of 
pressure.  These  are  depicted  in  Pig.  1  and  are  discussed  briefly 
below. 

The  first  regime  may  be  termed  the  Classical  Regime .  Here  the 
emphasis  is  on  quantities  which  vary  little  over  the  mean  free 
path  or  the  mean  time  of  intermolecular  collisions.  Moreover, 
molecular-surface  effects  are  quite  negligible  even  at  short  dis¬ 
tances  away  from  the  transmitter.  Indeed  the  theoretical  expres¬ 
sions  which  adequately  describe  sound  propagation  in  the  classical 
regime  do  not  contain  detailed  kinetic  parameters.  Nevertheless, 
the  use  of  the  kinetic  theory  to  derive  the  sound  propagation  in 
this  regime  serves  as  a  limiting  criterion  for  the  validity  of 
Eq.  (l.l)  and  may  also  provide  a  more  fundamental  understanding 
of  the  nature  of  sound. 

There  exists  a  range  of  parameters  which  we  designate  the  Frequency 
Relaxation  Regime .  Here  the  frequency  of  intermolecular  collisions, 
f  ,  is  comparable  to  or  smaller  than  the  frequency  f  of  the  sound 

V 

field  (f  >  f  ),  and  the  sound  field  is  observed  at  large  distances 

w*  w 


-2- 


Report  No.  1169 


Bolt  Beranek  and  Newman  Inc. 


from  the  transmitter  surface.  Then  sufficient  intemolecular 
encounters  have  taken  place  in  the  space  between  the  transmitter 
and  the  point  of  observation  t  0  ensure  that  the  observed  molecules 
possess  no  recollection  of  surface-molecular  interactions  so  that 
details  can  be  examined  entirely  on  the  basis  of  intermolecular 
collisions . 

In  the  Geometrical  Re laxatlon  Regime  the  molecules  are  assumed  to 
traverse  the  path  from  the  transmitter  to  the  receiver  control 
surface  essentially  without  undergoing  any  intermolecular  collisions. 

In  this  case,  one  sets  the  source  term  I _  of  Eq.  (1.1)  zero.  The 

solution  of  the  Boltzmann  equation  then  depends  primarily  on  the 
molecular-surface  interaction.  Consequently,  an  examination  of 
sound  propagation  in  the  geometrical  relaxation  regime  could  lead 
to  information  regarding  such  interactions. 

Our  chief  concern  is  with  the  third  regime  and  the  study  of  molecu¬ 
lar-surface  interactions.  The  research  reported  here  is  mainly 
exploratory  in  nature  in  both  its  experimental  and  theoretical 
aspects.  Nevertheless  some  of  the  findings  and  results  contribute 
in  some  measure  to  the  understanding  of  the  phenomenon  of  the 
propagation  of  collective  ordered  motion  in  a  rarefied  gas. 
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II.  DESCRIPTION  OP  THE  MODEL  FOR  DISPERSION  STUDIES 

To  avoid  complications  in  the  interpretations  of  both  the  experi 
mental  and  the  theoretical  results,  it  i3  expedient  to  impose 
certain  constraints  on  the  experimental  and  theoretical  models. 


1)  The  gas  molecules  only  possess  translational  degrees 
of  freedom  and  the  gas  is  composed  of  a  "single  com¬ 
ponent  . " 

2)  The  surfaces  of  the  transmitter  and  a  test  surface 
where  the  observations  of  the  sound  field  are  carried 
out  (e.g.,  hypothetical  control  surfaces  or  receivers) 
are  plane  and  parallel  to  each  other. 

3)  The  transmitter  surface  vibrates  uniformly  and  without 
distortion. 

4)  The  transmitter  is  large  enough  that  end  effects  are 
negligible . 

5)  Vibration  of  a  3ingle  frequency  is  superposed  on  the 
basic  motion  of  the  molecules. 

6)  The  amplitudes  which  define  the  sound  field  are  small 
compared  with  the  amplitudes  which  define  the  basic 
motion  of  the  molecules. 


We  shall  first  dwell  on  the  theoretical  aspect  of  the  phenomenon. 
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/ 

III.  DESCRIPTION  OP  DISPERSION  BY  THE  ATTENUATION 
AND  PHASE  (a  AND  0) 

Consideration  of  the  model  defined  above  permits  one  to  assume  a 
solution  to  the  equation  of  motion  of  the  form 

F(x,S,t)  -  F°(S)  +  h(x,$£)  ,  (3.1) 

w  ^ 

where  the  function  h  is  small  in  the  mean  (i.e.  the  ratio  of  a 
velocity  moment  in  h  to  one  in  F°  is  very  much  less  than  unity). 

In  Eq.  (3.1) 

x  *  position  of  a  control  surface  measured  in 
the  direction  normal  to  its  plane, 

£  =  molecular  velocity  vector, 

t  =  time  of  observation  at  the  control  surface, 

F°(£)  =  solution  to  the  equation  of  motion  in  the 
absence  of  sound. 

The  equation  of  motion  is  then  the  linearized  Poltzmann  equation 
which  is  stated  formally  in  Eq.  (l.l). 


The  sound  pressure  at  the  control  surface,  p(x,t)  is  the  normal 
flux  of  alternating  momentum  at  the  surface 


p(x,t) 


(3.2) 


where 


m  =  the  mass  of  each  molecule, 

£  =  a  unit  vector  in  the  x-direction. 
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The  dispersive  properties  of  the  sound  field  can  be  defined  by  the 
attenuation  per  unit  distance ,  a  and  the  wave  number  0,  for  waves 
at  angular  frequency  u>.  That  there  can  be  a  well  defined  frequency 
Is  evident  from  the  separability  of  Eq.  (1.1);  and  such  harmonic 
disturbances  can  be  Induced  by  Imposing  an  oscillating  boundary 
at  the  transmitter.  The  momentum  flux  (pressure)  at  any  space 
point  thus  is  of  the  form  p(x,t)  -  p(x)  exp(-iu>t).  However  no 
assumption  has  been  made  that  permits  one  to  assume  that  p(x)  is 
a  plane  wave  and  therefore  one  must  obtain  the  mathematical  defini¬ 
tion  of  a  and  0  from  an  analysis  of  the  measurement  process  for 
determining  attenuation  and  phase.  One  measures  the  pressure  at 
neighboring  space  points  In  the  vicinity  of  x,  obtains  the  difference 
and  normalizes  by  the  Interval  between  the  points  and  the  pressure 
itself.  Thus  one  obtains,  for  closely  neighboring  points  separated 
by  Ax, 

pl*  --<»-»>«  j$k  •  i«) 

Clearly  o  and  0  are,  in  general,  functions  of  x  and  therefore 
the  sound  pressure  partial  wave  of  frequency  u>  at  the  point  x  is 
to  be  expressed  in  terms  of  the  sound  pressure  at  x'  by  the  relation 

.  x 

p(x )  -  p(x')  exp  -  /  ds[a(s)  -  i0(s)]  .  (3.4) 

Jx' 

This  concept  of  a  generalized  attenuation  and  phase  which  are  not 
only  frequency  but  also  space  dependent  is  fundamental  to  the  subse¬ 
quent  discussion. 

Since  p(x,t)  is  a  moment  over  h(x,£,t)  we  must  first  derive  an 
expression  for  h(x,€,t)  in  terms  of  the  molecular  surface  interaction. 


-6- 


Report  No.  1169 


Bolt  Beranek  and  Newman  Inc 


IV.  DERIVATION  OP  A  GENERAL  SOLUTION 
A.  Formal  Solution 

Instead  of  attempting  a  direct  solution  of  Eq.  ( 1.1)  by  assuming 
a  suitable  form  for  the  source  term,  we  set  the  source  term  zero 
and  make  that  solution  subject  to  the  boundary  conditions  which 
the  source  term  repre  sents. 

l/e  denote  the  distribution  function  of  the  molecules  subsequent 
to  their  interaction  with  the  vibrating  transmitter  surface  by  G 
and  the  normalized  distribution  function  of  the  molecules  subse¬ 
quent  to  their  interaction  with  the  transmitter  surface  in  the 
absence  of  vibrational  motion  by  g_.  V/e  assume  that  the  incre¬ 
mental  velocity,  o,  and  density,  p°r},  which  is  imposed  on  the 
molecules  by  the  surface  vibration  are  linearly  superposed  on 
the  molecular  motion  in  the  absence  of  vibrations.  The  Darameter 
is  the  density  in  the  absence  of  vibrational  motion.  We  denote 
the  distribution  function  of  the  incremental  velocities  by  q.  G 
is  a  Joint  distribution  function  in  4'  and  o,  and  can  be  written 
in  product  form  because  J,'  and  o  are  independent.  Thus, 

mG_(x  =  0,^,o,t')  =  p°|l  +  T)(t'  )jq(o)g_(^’  -  o)U(^)  .  (4.1) 

Here  i *  =  molecular  velocity  vector, 

/  do  q(o)  =  1,  (4.2) 

j  d4g_($)  =1,  (4.3) 
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and 


0  ii  <  0 


f  'x 

x 

x  <•  1  ti  >  0 


(4.4) 


Ve  assume  that  1^1  /cm  and  i\  are  small  compared  with  unity  and  ex¬ 
pand  g_  In  Taylor's  series  retaining  only  the  first  order  terms  in 
0 


g-(i’  m  S)  m  «_(£')  -  ®  ^  e.(l' )  •  (4*5) 

Thus,  to  first  order  in  the  perturbation 

mCP  -  m  J 0_do  -  p°  |g_(£')  +  Tj(t')g_(^' ) 

’  S  *  3$T  g_<« • )}  W(«i)  »  (4.6) 


where 


o 


do  oq(o) 


and  r\ 


We  introduce  a  spatial  and  temporal  propagator  T  to  describe  the 
development  in  time  of  the  distribution  function  0^  as  x  goes 
from  zero  at  the  transmitter  to  its  value  at  the  control  surface. 
Formally  we  denote  this  development  by  the  equation 
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F(x,?, ;)  -  fit'  J  d?'  0;(x  -  0,  «',o,  t') 


T(x  0,£',  t ' |x,£,t)  • 


(4.7) 


The  determination  of  the  functional  form  of  T  la  equivalent  to 
solving  the  Boltzmann  equation  with  corresponding  constraints. 

This  usually  would  present  Just  as  many  difficulties  as  the  direct 
solution  of  Eq.  (1.1).  However,  in  the  geometrical  relaxation 
regime  the  propagator  assumes  the  simple  form 


T  *  6(t'  -  t  +  ■£-)  6(£'  -  £)  , 

*x  ~~  ^ 


(4.8) 


where  6  is  Dirac's  delta  function  and  6  .  668 

x  y  z 

Prom  Eqs.  (4.6),  (4.7)  and  (4.8)  one  obtains 

p  -  (p»  jg.(«)  +  n(t  -  £>  g  (?) 

-  S(t  -  •  &  «.(«)}  ®«x)  ’ 

and 

h  -  (P°A)  |qg_  “  Z  *  ^  g-J  U(^x)  • 


(4.9) 


Note  that  h  as  given  by  Eq.  (4.9)  satisfies  the  Boltzmann  equation 

with  both  I,  and  I  _  equal  to  zero  [see  Eq.  ( 1 • 1 ) ) •  The  function 
mm  ms 

h  satisfies  the  boundary  conditions  by  virtue  of  its  derivation 
from  Eq.  (4.7),  To  proceed  further  one  must  examine  these  boundary 
conditions  in  some  detail. 
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B,  Boundary  Conditions  (Molecular-Surface  Interaction 
at  the  Transmitter) 

Assume  that  the  distribution  function  of  the  incident  molecules 
on  the  surface  of  the  transmitter  is  (p^/in)g+(£  )U(-CX)  [p°  is 
independent  of  time  in  absence  of  vibrational  motion].  Consider 
the  transformation  from  coordinates  fixed  in  space  to  coordinates 
fixed  in  the  vibrating  surface  of  the  transmitter.  In  thi3  co¬ 
ordinate  system  the  distribution  function  for  the  incident  and 
emerging  molecules  is  given  by  (p£/m)g+(£  +  n)U(-Sx)  and(p  °/m) 

(1  +  hte./U  -  (0  -  £) ju(£x)  respectively.  The  velocity  u  is  the 
velocity  of  the  surface  of  the  transmitter  (note  that  u  =  &u). 

Now  suppose  that  the  surface  does  not  constitute  a  source  (or  a 
sink)  of  molecules.  One  can  then  equate  the  incident  and 
reflected  fluxes  and  obtain,  after  some  reduction,  the  conserva¬ 


tion  equations 

q  -  a(u/cm)  and  P°c_  -  P°c+,  (4.10) 

where 

a  “  cm[(v/2c+)  ‘  (y  "  1)/2c-^  *  (4.11) 

c+  -  J  di  *xg+  u(+^x)  '  (4.12) 

5X  -  YM-  ,  (4.13) 

y\L  m  j  dt'K(t' I  t)u(t' )  ,  (4.14) 

_00 


K( t * 1 1 )  ■  the  probability  that  a  molecule  that  strikes 
the  surface  at  time  t*  emerges  at  time  t. 
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Here  use  has  been  made  of  the  requirement  that  g+  (i)  — - •>  0 
as  |i|  — »  «  and  we  have  also  made  the  reasonable  assumption  that 
v  and  y  are  constants  In  time.  Note  that  v  and  y  depend  on  the 
type  of  molecular-surface  Interaction.  The  parameter  v  is  re¬ 
lated  primarily  to  the  "sitting  time"  of  the  molecules  on  the 
surface  of  the  transmitter.  This  surface  relaxation  time  can 
be  sensed  by  the  sound  field  when  the  average  sitting  time  is 
comparable  to  or  longer  than  the  period  of  the  sound  oscillation. 
The  parameter  y  is  related  primarily  to  the  efficiency  witn  which 
the  surface  impresses  its  vibrational  motion  on  to  the  molecules. 

The  exchange  of  the  normal  component  of  momentum  between  molecules 
and  surfaces  is  of  great  interest  to  aerodynamicists  and  may  be 
stated  in  terms  of  the  ratio,  PT,  of  reflected  to  incident  normal 
momentum  flux.  In  the  absence  of  surface  vibration  this  ratio  i3 

=  7c^c_  .  (4.15) 

In  terms  of  properties  of  the  sound  field  the  reflector  ratio  is 
P*  *  (y  -  1)  +  (ac3/2cmc_)  ,  (4.16) 


where 


c+  ■  J  d5  * 


and  a,  c+  and  y  are  defined  in  Eqs.  (4.11),  (4.12)  and  (4.13) 

respectively.  Equation  (4.16)  is  a  statement  that  acoustical 
measurements  could  be  used  to  ascertain  reflective  properties 
of  molecular-surface  systems. 
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C.  Classification  of  Surface  Interactions 

We  find  it  convenient  to  consider  the  molecular-surface  interaction 
in  terms  of  three  special  processes  by  which  the  general  case  may 
be  approximated.  Such  consideration  serves  to  elucidate  the  physi¬ 
cal  interpretation  of  the  quantities  v  and  y.  These  three  inter¬ 
action  processes  are  outlined  below: 

1.  Specular  Reflection 

Each  molecule  is  assumed  to  interact  with  the  surface  instantly 
and  in  such  a  way  that  the  sign  of  its  component  of  momentum  in 
the  x-direction  is  reversed.  This  interaction  process  conserves 
energy  and  momentum.  It  is  easy  to  show  that  in  this  case  (viz. 

Eq.  (4.1)  et  seq. ) 

y  =  2  and  v  =  1 

2.  Diffuse  Elastic  Reflection 

Each  molecule  is  assumed  to  interact  instantly  with  the  surface 
3uch  that  the  incident  energy  is  conserved  but  not  the  momentum. 

It  is  further  assumed  that  the  scattering  obeys  the  cosine  law 
of  reflection.  In  this  case  one  can  3how  that 

y  =  ^  and  v  -  1 

3.  Diffuse  Inelastic  Reflection 

The  molecules  are  considered  adsorbed  by  the  surface  and  to  come 
to  thermal  equilibrium  with  it  before  boiling  off  in  a  random 
fashion.  That  is,  it  is  assumed  that  neither  the  momentum  nor 
the  energy  are  conserved,  and  that  the  distribution  of  the 
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molecules  subsequent  to  the  interaction  is  consistent  with  the 
assumption  that  the  gas  is  in  thermal  equilibrium  with  .he  surface. 

In  this  case  there  are  two  regimes.  If  the  function  K(t'|t) 
is  such  that  most  of  the  contribution  to  the  integral  in 
Eq.  (4.14)  comes  from  a  region  where  t  -  t'  <  then 

v  =  1;  in  the  other  extreme  v  =  0.  For  diffuse  inelastic  re¬ 
flection  7  =  1. 

When  more  complicated  molecular-surface  interactions  can  be 
approximated  by  such  elementary  processes,  the  contribution 
from  each  basic  process  is  computed  separately  and  added  to 
give  the  total  effect. 
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V.  DETERMINATION  OF  THE  TRANSMISSION  OF  SOUND  IN 
TOE  GEOMETRICAL  RELAXATION  REGIME 

The  equation  of  sound  propagation  in  the  geometrical  relaxation 
regime  can  be  determined  using  the  semi-phenomenological  parameters 
introduced  in  the  previous  section.  Substituting  Eqs.  (4.9)  and 
(4.10)  in  Eq.  (3.2)  and  maldng  straightforward  manipulations,  we 
obtain 


00 

p(s)  -  P>0°m  Jo  dr  1,2 

x  SgXl 

■jag_(r)  -  y  J  exp(is/r)  , 

(5.1) 

where 

3  *  J*/cra  > 

(5.2) 

r  *  Vcm  • 

(5.3) 

Note  that  Eq.  (5.1)  is  dependent  on  the  molecular-surface  inter¬ 
actions  in  a  complex  way.  There  are  cases,  however,  for  which  it 
is  possible  to  assume  the  analytical  form  of  g+  and  g  ;  for  such 
cases,  v  and  y  could  be  examined  by  acoustical  means.  Conversely, 
if  v  and  Y  could  be  estimated  by  an  Independent  method,  the  distri¬ 
bution  function,  g  ,  could  be  studied  by  acoustical  means. 


Equation  (5.1)  is  suggestive  of  the  work  of  other  investigations. 

It  differs  from  the  equation  derived  by  Meyer  and  Sessler^but 
for  s  »  1  and  for  a  Maxwellian  distribution  it  yields  the  same 
dispersion  equations.  Our  results  agree  with  an  equation  derived 
by  Kahrv^  when  the  following  restrictions  are  imposed  on  Eq.  (5.1): 
the  distribution  functions  arc  Maxwellian;  v  *  0  and  y  =  1  (i.e., 
diffuse  inelastic  reflection  and  t  -  t'  of  Eq.  (4.14)  is,  on  the 
average,  very  much  larger  than  7r/2to). 
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VI.  PHYSICAL  INTERPRETATION  OF  SOUND  DISPERSION 
IN  THE  GEOMETRICAL  RELAXATION  REGIME 

From  Eq.  (5.1)  It  is  evident  that  the  expression  for  the  sound 
pressure  has  the  general  form 

oc 

p(s)  =  N  /  dr  A(r)  exp  (is/r)  ,  (6.1) 

Jo 

In  which  N  is  independent  of  s.  It  is  through  a  careful  analysis 
of  the  significance  of  the  various  contributions  to  the  integral 
in  this  expression  that  we  uncover  the  mechanism  whereby  the 
sound  wave  suffers  dispersion  in  propagation. 

The  integral  has  an  appearance  very  similar  to  the  Fourier  integral 
and  as  such  can  be  thought  of  as  being  a  superposition  of  weighted 
plane  waves.  A(r)  would  be  the  weighting  function  and  the  phase 
term  exp  is/r  a  plane  wave  of  wave  number  ^  at  dimensionless  posi¬ 
tion  s.  The  weighting  function  A(r)  is  the  product  of  the  molecular 
velocity  distribution  function  and  a  factor  which  brings  the  tenns 
N  A(r)  dr  to  the  dimension  of  momentum  flux,  l.e.,  pressure.  Thus 
the  partial  sound  pressure  (the  sound  pressure  contribution  from 
an  infinitesimal  velocity  range)  i3  given  by 

dp(s)  *  NA(r)exp(is/r)dr  .  (6.2) 

From  the  above  discussion  and  Eq.  (6.2)  one  can  conclude  that 
were  there  but  one  velocity  there  would  be  no  dispersion  and 
the  wave  number  would  be  inversely  proportional  to  that  one 
velocity. 
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Extrapolation  leads  us  to  state  that  for  a  highly  peaked  flux 
distribution  there  would  be  little  attenuation  and  dispersion. 

In  nature,  however,  the  fluid  molecules  have  a  rather  broad 
velocity  distribution.  For  equilibrium  and  near  equilibrium 
the  weighting  function  A(r)  has  the  general  form  depicted  in 
Fig.  2a.  There  is  a  maximum  value  of  A  when  r  is  of  the  order 
of  unity.  Consider  now  the  influence  that  such  a  broad  weighting 
function  has  on  a  superposition  of  the  partial  waves  of  Eq.  (6.2). 
In  the  superposition  the  partial  waves  can  interfere  destructively 
or  constructively.  The  kind  of  interference  contributed  in  the 
various  ranges  of  the  variable  r  in  the  integration  can  be  dis¬ 
cussed  by  considering  s/r  in  the  ranges  3/r  >  1  and  s/r  <  1. 

The  behavior  of  the  phase  term  as  a  function  of  the  variables  r 
and  s  is  depicted  in  Figs.  2b,  c  and  d.  For  values  of  s/r  large 
compared  with  unity  the  phase  term  is  an  oscillatory  function  of 
the  variable  r;  the  slower  molecules  (s  >  r)  arrive  at  the  re¬ 
ceiver  with  widely  varying  phases.  The  partial  sound  pressures 
associated  with  these  molecules  interfere  with  one  another 
destructively  (phase  mixing).  Thus  the  slower  molecules  do  not 
contribute  to  the  sound  pressure.  The  faster  molecules  (3  <  r) 
arrive  at  the  receiver  with  essentially  the  same  phase.  This 
preservation  of  coherence  by  the  faster  molecules  constitutes 
the  sound  field  at  the  receiver.  The  inability  of  the  slower 
molecules  to  transmit  sound  is  directly  responsible  for  the  dis¬ 
persive  properties  of  sound.  Thus,  what  constitutes  slow  or  fast 
molecules  is  determined  by  whether  the  molecules  possess  velocity 
components,  r  smaller  or  greater  than  s,  respectively.  To  account 
for  the  dispersive  properties  of  sound  in  more  detail  it  is  there¬ 
fore  convenient  to  consider  the  problem  in  three  ranges  of  the 
parameter  s.  They  are  s  «  1,  s  »  1  and  s  *  1  (see  Fig.  2). 
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In  order  to  introduce  some  quantitative  data  in  the  following 
discussion  reference  is  made  to  the  example  of  a  Maxwellian 
distribution.  That  is,  we  assume  that  both  g+  and  g_  are 
Maxwellian  distributions.  Then  from  the  foregoing  equations  we 
obtain 

N  =  p  u  c  /fir, 
o  o  m' 

A(r)  =  r2(a  +  2yr)  exp  (-r2)  , 

where 

p  =  U0  exp  (-i“>t). 


and 

a  =  Vtt(v  +  1  -  Y)  . 

A.  Dispersion  in  the  Range  s  «  1;  Figs.  2a  and  2b 

In  this  range  only  the  very  slow  molecules  do  not  contribute 
to  the  sound  transmission.  Since  their  associated  weighting 
functions  are  small  they  do  not  impair  greatly  the  total  transfer 
of  coherence  from  the  transmitter  to  the  receiver.  Thus  one  would 
expect  the  attenuation  parameter  a  to  be  small  and  to  increase 
monotonically  with  increase  in  s.  The  phase  parameter  p  (wave 
number),  which  is  inversely  proportional  to  the  speed  of  sound, 
is  not  greatly  affected  in  this  range  by  variation  in  s.  As  s 
increases  a  corresponding,  but  slight,  decrease  in  $  is  expected 
since  it  is  left  to  the  faster  molecules  to  transmit  the  sound 
(coherence) . 


(6.3) 

(6.4) 
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We  cite  the  results  which  are  computed  by  solving  Eq.  (3.4)  using 
Eq.  (6.1)  and  the  Maxwellian  quantities  given  in  Eqs.  (6.3)  and 
(6.4) 

(a  -  i£)/P0  -  saQ[l  +  a^2  +  ...) 

-  it>o [  1  +  b^32  +  ...]  1  (6.3) 

where 

aQ  ■  (^)  '  [{sJir  +  Uy)(st/v  +  2y) 

-  2(a  +  Wy)2]/(Vira  +  hy)2  ,  (6.6) 

m  V2 

bQ  «=  (- j )  (a  +  VrTry)/(Vira  +  4y)  .  (6.7) 


*  the  phase  parameter  for  an  ideal  gas  under 
standard  condition  of  temperature  and  pressure, 
al***bl  “  constants  which  depend  on  v  and  y. 

Note  that  over  reasonable  limits  on  v  and  y  (0  <  v  <  l  and 
1  <  y  <  2]  the  value  of  bQ  is  near  unity. 

The  quantitative  resu?it,  Eq.  (6.5)>  bears  out  the  qualitative 
discussion. 
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B.  Dispersion  In  the  Range  s  »  1;  Figs.  2a  and  2d 

In  this  range  only  the  very  fast  molecules  can  participate  in 
transmitting  sound.  Thus  the  phase  parameter  £3  is  small  compared 
with  that  in  the  range  s  «  1.  However,  its  rate  of  decrease 
with  increase  in  s  is  expected  to  be  proportionately  slow.  The 
reason  is  that  the  weighting  function  is  a  monotonically  decreas¬ 
ing  function  of  r;  thus  the  molecules  which  are  concentrated  in 
the  vicinity  of  r  -  3  are  those  which  most  influence  the  trans¬ 
mission.  However,  the  cessation  of  the  oscillation  occupies  a 
relatively  wide  range  in  the  parameter  s/r.  Both  these  effects 
combine  to  resist  large  changes  in  3  as  a  function  of  s.  The 
attenuation  parameter  a[N.B.  a  i3  a  local  attenuation  variable, 
see  Eq.  (3.3)]  follows  similar  changes  with  s  as  does  the  phase 
parameter,  that  is,  its  value  decreases  at  a  proportionate^  slow 
rate  with  increase  in  s. 

Again  we  make  use  of  the  Maxwellian  distribution  to  obtain 

quantitative  results  for  this  range.  Making  use  of  the  method 

4 

of  steepest  descent-'  one  obtains 

(a  -  ie)/P0  =  (5!) 1/2 (5)' 1/3  {l  -  W3 f  ,  (6.8) 

which  is  in  agreement  with  the  preceding  qualitative  discussion. 

The  dispersion  in  this  range  is,  to  first  order  of  approximation, 
devoid  of  molecular-3urface  interaction  parameters.  The  results, 

Eq.  (6.8),  are  given  in  graphical  form  in  Pig.  3. 

C.  Dispersion  in  the  Range  s  -  l;  Flgc.  2a  and  2c 

In  the  range  r  ■*  1  the  weighting  function  is  relatively  large  and  it 
undergoes  large  variation  as  a  function  of  r.  Thus  one  would 
expect  large  variation  to  occur  in  both  the  dispersion  parameters. 
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The  attenuation  parameter  Increases  as  s  Increases  to  unity.  As  s 
Is  Increased  further,  the  attenuation  reaches  a  maximum  value  and 
then  decreases  to  an  asymptote  for  s  »  1.  Ihe  phase  parameter 
decreases  significantly  with  Increase  In  s  only  when  s  becomes 
greater  than  unity.  For  only  then  Is  the  weighted  mean  velocity 
component  of  the  faster  molecules  In  the  direction  of  propagation 
Increased. 

In  this  range  we  were  not  able  to  obtain  analytical  expressions 
for  p(s)  or  its  derivative  with  respect  to  s,  therefore,  we  per¬ 
formed  numerical  computations.  The  results  are  plotted  in  Fig.  3. 
Observe  that  the  dependence  of  the  dispersion  on  the  molecular- 
surface  interaction  parameters  is  present  but  not  pronounced. 

This  Is  unfortunate;  it  indicates  that  considerable  accuracy  in 
the  measurements  of  the  dispersion  of  sound  is  required  to  obtain 
results  that  can  distinguish  between  the  various  processes  of 
surface-molecular  interactions . 
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VII.  EXPERIMENTAL  WORKING  RANGE 

In  order  to  extend  the  range  of  the  work  which  has  been  previously 
conducted  in  this  field, as  well  as  to  examine  new  applications 
where  sound  nropagation  may  be  useful  for  studies  of  rarefied  gas 
dynamics,  we  chose  to  conduct  our  experiments  at  lower  frequencies* 

(2  to  12  kc)  and  at  lower  pressures*  (as  low  as  5  x  10"^  mm  Hg). 

The  lower  frequency  range  has  the  advantage  that  sound  measurements 
can  be  made  in  a  truly  rarefied  gas  over  larger  separations. 
Measurements  over  longer  distances  have  the  advantages  that  higher 
accuracies  can  be  achieved  and  that  auxiliary  equipment  can  be 
chosen  without  the  pains  of  miniaturization.  The  challenge  to 
this  approach  is  whether  one  can  produce  high  enough  signals  and 
sensitive  enough  receivers  to  achieve  measurements  with  low  sig- 
nal-to-noise  at  the  lower  pressures. 


‘l^atively 


Greenspan-^  conducted  his  measurements  at  11  Me  and  at  r$ 
high  pressures,  higher  than  1  mm  Hg.  Meyer  and  3ecsler£/con- 
ducted  their  experiments  at  100  and  200  kc  and  at  pressures  in 
excess  of  2  x  10“^  mm  Hg. 
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VIII.  THE  EXPERIMENTAL  SYSTEM  FOR  STUDIES  OF  DISPERSION 

In  this  section,  we  describe  briefly  the  experimental  system 
which  we  used  to  study  the  dispersion  of  sound  in  rarefied  air. 

The  system  is  designed  to  operate  in  pressure  ranges  where  the 
mean  free  path  of  intermolecular  collisions  is  longer  than  the 
spacing  between  the  transmitter  and  receiver. 

A.  The  Vacuum  System 

For  good  accuracy  and  flexibility,  it  is  desirable  that  the 
spacing  between  the  transmitter  and  receiver  be  varied  over 
several  centimeters.  Separations  of  the  order  of  1  cm  or  more 

_  'i 

require  working  pressures  of  the  order  of  10  mm  Hg  or  lower. 

(In  air  the  mean  free  path  at  10"  ^  mm  Hg  and  300°  K  is  about 
5.5  cm.)  Since  mechanical  pumps  cannot  achieve  these  low  pres¬ 
sures,  it  is  necessary  to  use  a  diffusion  pump. 

The  pumping  system,  which  is  of  conventional  design,  is  shown  in 
Fig.  4.  The  limiting  pressure  obtained  in  the  test  section  with 
this  system  wa3  about  3  x  10" 0  mm  Hg. 

The  transmitter-receiver  and  the  bell  jar-base  plate  systems  are 
shown  in  Fig.  5.  The  sound  source  (transmitter)  is  mounted  on 
a  movable  vacuum  tight  rod  which  can  be  actuated  externally. 

With  this  arrangement,  the  distance  between  the  source  and  receiv¬ 
ing  microphone  can  be  varied  continuously  over  a  r  .ge  of  12  cm 
while  maintaining  a  constant  pressure  in  the  test  section. 

Two  vacuum  gauges  are  incorporated  in  the  system.  A  thermo- 
couple  gauge  is  used  in  the  range  from  10  mm  Hg  to  5  x  10  J  mm 
Hg  and  an  ionization  gauge  covers  the  range  from  5  x  10"**  mm  Hg  to 
10"6  mm  Hg.  We  also  had  available  an  alphatron  gauge  with  a  range 
from  10  mm  Hg  to  10  J  mm  Hg. 
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E.  The  Electrical  and  Mechanical  System 

The  microphone  is  of  the  electrostatic  type  (Bruel  and  KJaer 
Model  4132).  A  1-inch  microphone  is  suspended  vertically  above 
the  sound  source  in  a  fixed  location.  It  is  mounted  in  such  a 
way  that  the  electrical  connections  are  through  the  top  of  the 
bell  Jar.  This  arrangement  is  chosen  in  order  to  reduce  the 
background  signals  due  to  electrical  stray  currents.  The  polari¬ 
zation  voltage  of  the  microphone  is  set  so  that  no  discharge 
occurs  between  membrane  and  backing  plate  at  the  operating  pres¬ 
sure  . 

The  sound  source  is  an  electrostatic  transducer  with  solid 
dielectric.  The  vibrating  diaphragm  is  light  enough  so  that 
mechanical  cross  talk  between  it  and  the  microphone  is  essentially 
eliminated.  The  vibrating  surface  extends  to  about  a  5.0  inch 
diameter.  The  output  of  the  sound  source  has  a  maximum  of  130  db 
at  atmospheric  pressure. 

T  .e  electric*  1  diagram  of  the  experimental  apparatus  is  shown  in 
Fig.  6. 

The  separation  between  microphone  and  transmitter  can  be  varied 
automatically,  continuously  and  slowly  over  a  range  of  about  12  cm. 
The  initial  separation  may  be  set  as  small  as  0.05  cm.  The  separa¬ 
tion  x  as  a  function  of  time  i3  given  by 

x  =  b  +  d [ 1  -  cos  (t/r)]  ,  (8.1) 

where  b  is  the  initial  separation  (which  can  be  varied  in¬ 
dependently),  d  =  12.1  cm  and  r  is  set  at  about  90  sec¬ 
onds  (t  may  also  be  controlled  independently).  This 
particular  rate  of  change  of  stroke  was  decided  upon  because 
it  is  simple  to  Implement  and  also  enables  one  to 
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expand  the  Initial  region  (t  «  r)  where  the  attenuation  and  phase 
parameters  show  definite  deviation  from  the  theory  in  which  mul¬ 
tiple  reflections  are  neglected.  The  disadvantage  of  this  varia¬ 
tion  is  that  the  records  have  to  be  converted  to  a  linear  scale. 

We  have  therefore  installed  mechanical  and  electrical  linear 
scale  converters  which  can  be  "switched  in"  when  records  over 
large  separations  are  desired.  With  these  converters  operating, 
the  abscissa  of  the  graphic  record  is  proportional  to  the  distance 
travelled  by  the  transmitter.  With  these  devices,  it  is  estimated 
chat  the  spacing  can  be  read  to  better  than  0.1  cm  over  the  entire 
range  of  separation. 

The  sound  pressure  level  13  recorded  on  a  graphic  level  recorder 
and  can  be  read  tc  within  0.5  db  on  a  30  db  potentiometer.  When 
a  limited  range  of  the  sound  pressure  level  is  of  interest  a 
10  db  potentiometer  is  used  and  the  accuracy  thereby  improved. 

Tke  results  which  are  reported  here  were  obtained  with  a  50  db 
potentiometer  only.  Typical  records  are  reproduced  in  Fig.  7. 

The  ldentatlcn  in  Fig.,  7a  at  about  1  cm  is  due  to  an  approxi¬ 
mately  10  db  change  in  gain  which  is  sometimes  necessary  to  keep 
the  entire  record  on  a  50  db  potentiometer. 

The  phase  is  measured  with  a  phasometer  which  indicates  the  phase 
difference  between  a  direct  signal  from  the  driving  oscillator 
and  the  signal  received  at  the  microphone.  The  output  of  the 
phasometer  is  fed  into  the  ordinate  of  a  graphic  recorder.  We 
estimate  that  in  most  cases  the  phase  could  be  read  to  within 
10°  over  a  phase  change  of  360°.  A  typical  record  is  shown  in 
Fig.  8.  The  Jump3  in  the  curves  (at  2  cm  and  around  10  cm) 
occur  when  one  full  meter  scale  reading  (360°)  is  completed. 
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The  useful  frequency  range  of  the  present  apparatus  extends  from 
5  kc  to  about  12  kc.  The  upper  limit  is  set  by  the  sensitivity  of 
the  microphone  which  drops  off  considerably  beyond  its  resonance 
at  about  10  kc.  The  lower  limit  i3  set  by  the  available  stroke 
of  12  cm.  Lower  frequencies  are  useful  for  studies  of  multiple 
reflections  only. 
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IX.  EXPERIMENTAL  RESULTS 

All  of  the  many  records  which  were  taken  (e.g.,  Pigs.  7,  8,  and  9) 
show  that  there  are  two  distinct  regimes  in  the  behavior  of  the 
amplitude  and  phase  of  the  transmitted  sound.  There  13  the 
"initial  region"  (up  to .approximately  1  cm)  where  the  signal 
decreases  very  rapidly  and  the  sign  of  the  phase  is  opposite 
to  that  predicted  by  the  theory.  It  is  only  at  a  separation  of 
more  than  1  cm  thft  the  attenuation  and  the  phase  of  the  signal 
follow,  at  least  qualitatively,  the  theory  outlined  above. 

The  region,  where  multiple  reflections  can  be  ruled  out  (x  »  1  cm), 
we  term  the  "far  field.  The  region  between  these  two  extremes 
we  refer  to  as  the  "intemediate  region."  Roughly,  the  initial 
region,  the  intermediate  region  and  the  far  field  correspond  to 
regions  which  are  defined  by  s  «  1,  s  1  and  3  »  1  respectively, 
where  s  is  the  normalized  separation.  We  shall  now  consider  each 
region  separately. 

A.  The  Initial  Region 

We  found  that  the  initial  decay  with  increasing  separation  is 

_  k 

independent  of  the  pressure  down  to  5  x  10  J  mm  Hg  for  a  given 
frequency,  and  up  to  a  separation  of  about  1  cm.  A  slight  varia¬ 
tion  of  the  initial  decay  with  variation  in  frequency  is  observed. 

We  conclude  that  for  small  separations  (<  1  cm)  and  at  a  given 
frequency  the  decay  with  distance  is  independent  of  the  pressure 
(this  was  observed  over  pressures  ranging  from  5  x  10”"  mm  Hg  to 
10  mm  Hg)  in  the  frequency  range  from  2  to  8  kc .  To  establish 
tnat  the  initial  decay  is  not  a  result  of  a  possible  variation  in 
amplitude  and  phase  over  the  surface  of  the  transmitter,  the 
phenomenon  was  examined  at  various  sections  of  the  transmitter. 

The  results  of  these  experiments  show  no  variation  in  the  detailed 
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structure  of  the  Initial  decay  curve.  Furthermore,  we  replaced 
the  electrostatic  transmitter  by  a  5-inch  plane  piston  radiator 
(first  resonance  at  about  15  kc )  which  was  driven  by  aui  electro¬ 
mechanical  shaker.  We  compared  the  initial  decay  curves  obtained 
at  various  pressures  and  found  them  to  be  the  same  within  the 
accuracy  of  the  measurements.  The  piston  radiator  was  only  used 
in  the  lower  frequency  range  (<  3  kc )  due  to  its  poor  efficiency 
at  higher  frequencies. 

Using  our  theory,  we  are  able  to  partially  account  for  the  be¬ 
havior  of  the  sound  in  this  region.  That  the  slope  of  the  sound 
pressure  level  curve  is  independent  of  pressure  at  small  separa¬ 
tions  may  be  explained  when  one  considers  that  for  small  separa¬ 
tions,  the  attenuation  and  phase  in  the  geometrical  relaxation 
regime  resemble  the  corresponding  parameters  in  the  classical 
regime.  This  resemblance  is  not  only  qualitative  but  also 
quantitative  (see  Section  VI).  Since  the  attenuation  and  phase 
define  the  transmission  of  sound  and  since  we  see  no  reason  why 
the  presence  of  multiple  reflections  should  alter  this  corres¬ 
pondence,  no  transition  region  is  to  be  expected  as  one  proceeds 
from  the  classical  regime  to  the  geometrical  relaxation  regime. 

B.  The  Intermediate  Region 

In  this  region,  commencing  first  at  higher  frequencies  (c.f. 

Fig*  9)*  the  multiple  reflections  lade  and  the  response  is  pre¬ 
dominately  the  response  of  a  single  reflection  at  each  surface. 
Figures  7  and  9  (which  are  representative  figures)  show  that  the 
attenuation  is  relatively  small  in  this  region.  However,  it 
increases  somewhat  with  increasing  separation.  Comparison  with 
Fig.  3  indicates  that  this  behavior  is  in  agreement  with  theory. 
The  phase  also  assumes  a  behavior  which  is  more  commensurate  with 
the  single  reflection  theory. 
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C.  The  Far  Field  Region 

In  thia  region,  (x  >  2  cm)  we  can  test  the  theory  which  we  derived 
above.  The  experimental  results  are  plotted  in  Pigs.  10  and  11. 

Aa  abaclaaa  in  theae  figures,  we  use  the  quantity  R  *  (2/3ir)(f  /f)|/ 
It  is  related  to  the  gas  pressure  p  through  f  ,  the  mean  frequency 

w  C 

of  intermolecular  collisions.  A  corresponding  parameter  can  be 
defined  for  the  geometrical  relaxation  regime,  ^  -  2/3 ir  c^fxf) 

■  4/(38).  Then  cm/x  is  the  mean  frequency  of  molecular-surface 
collisions.  Note  that  in  the  geometrical  relaxation  regime  the 
attenuation  and  phase  are  independent  of  the  collision  frequency 
f  and  hence  of  the  ambient  pressure. 

v 

Curves  predicted  by  the  classical  Burnett  theory  and  the  values 
predicted  by  our  theory  (in  the  geometrical  relaxation  regime) 
are  also  plotted.  The  computations  are  based  on  the  diffuse 
elastic  process  of  molecular- surface  interaction.  There  is  an 
excellent  agreement  between  theory  and  experiment  in  the  geometri¬ 
cal  relaxation  regime.  The  experimental  data  which  are  analyzed 
for  values  of  the  dimensionless  separation  parameter  s  =  5.22  and 
s  •  10.44  seem  to  favor  the  diffuse  elastic  process  of  molecular- 
surface  interactions  (see  Fig.  3).  However,  more  refined  experi¬ 
ments  and  analyses  of  the  data  would  be  required  before  a  definite 
conclusion  on  the  type  of  surface  interaction  is  Justified. 

The  transition  between  the  frequency  relaxation  regime  and  the 
geometrical  relaxation  regime  is  of  special  interest.  If  one 
accepted  the  Burnett  theory  as  fairly  representative  of  the  situa¬ 
tion  in  the  classical,  and  in  the  frequency  relaxation  regime 
(even  in  the  range  of  pressures  where  Au>/cm  <  2)  then  one  might 
conclude  (from  Pig3.  10  and  11)  that  the  transition  region  from 
the  frequency  relaxation  regime  to  the  geometrical  relaxation 
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regime  extends  over  a  pressure  range  of  about  half  a  decade.  This 
Is  particularly  evident  in  Fig.  10  where  this  transition  region 
extends  into  the  range  of  pressure  for  which  the  attenuation  has 
its  maximum  value.  The  experimental  data  here  have  values  which 
fall  somewhere  between  the  Burnett  predictions  and  those  from  the 
theory  which  applies  to  the  geometrical  relaxation  regime  (c.f. 
Fig.  3). 

The  predictions  of  the  present  theory  are  based  on  the  assumption 
that  the  velocity  distribution  of  the  molecules  subsequent  to 
their  interaction  with  the  surface  is  Maxwellian.  We  have  not 
ascertained  the  magnitude  of  the  deviation  from  this  distribution 
that  is  required  to  significantly  influence  the  results,  however 
there  is  consistency  between  this  assumption  and  the  experimental 
data. 
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X.  SOUND  PROPAGATION  IN  RAREFIED  OASES  —  WAVE  CONCEPT 

In  thl3  section  a  wave  theory  approach  to  the  interpretation  of 
sound  propagation  in  rarefied  gas  is  briefly  expanded.  The  pur¬ 
pose  of  this  analy3i3  is  to  construct  an  alternative  analysis  to 
the  one  of  the  particle  concept  that  was  given  above.  This  is 
brought  out  not  only  to  reassess  the  results  obtained  by  the 
particle  concept  but  also  in  the  hope  that  some  aspects  of  the 
problem  that  cannot  easily  be  interpreted  in  terns  of  the  one 
may  be  more  readily  interpreted  by  the  other.  In  particular, 
the  phase  reversal  that  occurs  in  the  initial  region  is  explained 
more  readily  by  the  wave  concept. 

A.  The  General  Idea 

First  consider  only  those  molecules  that  have  the  3ame  thermal 
222  1/2 

speed  c  -  (i  +  i  +  £  )  *  and  assume  specular  reflection. 

a  y  z 

Then  the  speeds  before  and  after  interaction  with  the  vibrating 
surface  are  equal  (except  for  the  small  ac  component),  and  the 
molecules  leaving  a  surface  element  are  equally  distributed  over 
all  angles.  The  molecules  which  emerge  from  a  surface  element 
have  properties  which  are  the  same  as  that  of  a  wave  emanating 
from  a  surface  element.  They  have  a  given  3peed  of  propagation 
and  they  are  equally  distributed  over  all  angles.  Therefore, 
the  additional  average  momentum  that  is  given  to  the  molecules 
by  a  vibrating  surface  element  can  be  considered  as  a  wave 
propagating  out  from  this  surface  element.**  Hius  by  treating 


We  impose  an  isotropic  velocity  distribution. 

This  13  not  the  usual  particle-wave  analogy  where  the  particle 
itself  is  represented  by  a  wave  whose  wavelength  is  given  by 
the  total  momentum  of  the  particle.  In  our  case,  only  the 
additional  momentum  represents  a  wave. 
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the  molecules  carrying  coherent  momentum  from  the  surface  by  the 
technique  of  wave  propagation  analysis  one  makes  use  of  the 
elements  of  classical  acoustics  to  Interpret  sound  propagation  in 
a  rarefied  atmosphere.  There  is  also  the  convenience  that  the 
angular  distribution  is  averaged  out  initially. 

P .  Radiation  from  a  Moving  Surface 

Still  considering  only  molecules  with  the  same  thermal  speed  c 
take  the  velocity  of  the  moving  surface  element  dS  to  be  u  and 
define  a  wave  number  k  =  ^o/c  (^  =  angular  frequency  of  the  driver). 
The  elemental  sound  nressure  dp  that  is  generated  by  the  surface 
element  do  at  a  distance  h  is  theiW 

dp  .  .ipuJ  23^  ^  .  (10.1) 

(;  =  density  of  tne  gas.) 

The  result  for  the  total  vibrating  surface  S  is  obtained  by 
integrating  over  dS.  To  obtain  the  actual  sound  pressure  in  a 
rarefied  gas  all  v  ilues  of  c  must  be  considered.  This  is  ac- 
connlished  by  multinlying  Eq.  (10.1)  by  the  distribution  function 
f(c)  for  the  thermal  speed  c  md  integrating  over  dc .  Thus 


00 


P  = 


-up 

'  JOi.) 


.  F(e)  (led:; 


For  a  gas  at  c<-iiliL~  urn 


F(c)  =  ^  exnf - ( c  cjd) 

V  7T  C 

m 


(10.2) 


(10. j) 
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Therefore 


i-y. 


P  -  -7-  u  /  exp[-(c/o  )  ] 

/r  Jo  cf  m 

m 


{fs  riaf  *%(lkh)  dS)  de  • 


(10.4) 


Thu3  for  any  configuration  for  which  the  radiation  behavior  (the 
integral  over  S)  is  known  the  propagation  properties  can  be  found 
by  integrating  over  the  speeds  of  propagation.  This  is  especially 
useful  if  one  wishes  to  calculate  the  sound  propagation  from  a 
complicated  transmitter  system  in  a  rarefied  gas.  To  apply 
Eq.  (10.4)  one  calculates  the  wave  propagation  from  the  geometry, 
multiplies  it  by  the  velocity  distribution  function  and  integrates 
over  all  possible  values  of  c. 

C .  Radiation  from  an  Infinite  Disk 

For  radiation  from  an  infinite  disk  the  integration  over  S  is 
readily  performed.  The  integral  in  the  brackets  yields  the 
expected  plane  wave.  One  obtains 


=  pc  exp  iloc  =  pc  exP(lu*/c) 


(10.5) 


in  which  x  denotes  the  normal  distances  from 
infinite  di3k  to  receiver. 
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By  Inserting  Eq.  (10.5)  into  Eq.  (10.4)  one  obtains  for  the  pressure 


P 


exp[-r 


2 


+  i(s/r)]  dr 


(10.6) 


where 


r  *  c/e  and  s  -  u&/c  . 
m  'm 

This  is  exactly  Eq.  (5.1)  for  specular  reflection,  (v  ■  1,  y  -  2.) 

Since  p  and  cm  are  known,  a  measurement  of  u  together  with  a  known 
solution  of  the  integral  (i.e.,  for  s  >  3)  provides  a  method  of 
checking  numerical  constants  in  Eq.  (10.6).  Such  a  measurement, 
if  its  accuracy  is  better  than  1  db,  would  indicate  whether 
specular  or  diffuse  reflection  is  dominant.  This  accuracy  could 
be  attained  only  by  refinement  of  the  above  experimental  design. 

D.  Nonspecular  Reflection  from  the  Transmitter 

The  above  calculations  included  only  the  case  of  specular  reflec¬ 
tion,  Other  reflection  properties  could  also  be  considered  by 
using  the  "wave  concept."  However,  to  do  this  one  would  have  to 
know  the  radiation  properties  of  materials  with  finite  impedance 
(how  much  sound  is  radiated  from  a  vibrating  and  absorbing  sur¬ 
face?).  Unfortunately,  no  complete  solution  of  this  problem  is 
available.  We  assume  that  impedance  of  a  moving  disk  alters 
mainly  the  amplitude  of  the  radiated  sound;  other  effects  seem  to 
be  less  important.  (This  seems  to  hold  for  diffuse  elastic  and 
diffuse  inelastic  reflections.)  Thus,  we  continue  to  confine  the 
discussion  to  specular  reflection  only  in  using  the  wave  concept. 
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E.  A  Finite  Microphone  In  Front  of  an  Infinite  Transmitter 

(Sound  Waves  under  Standard  Conditions) 

Consider  the  transmitter-microphone  system  depicted  in  Pig.  12. 

The  classical  sound  pressure  pQ  as  measured  by  the  microphone 
can  be  expressed  by 

p  .  upc  t1  ♦  g)  e*p  .  (10.7) 

3  c  1  ■  R  exp  2ikx 

(u  *  velocity  of  the  transmitter,  k  =  wave  number,  R  *  reflection 
coefficient. ) 

Equation  (10.7)  is  obtained  by  adding  up  all  the  reflected  waves. 

Note  that  for  R  «  1,  one  obtains  the  usual  plane  wave  solution. 

o*r 

R  is  a  complicated  function  of  x,  d  and  X  =  ,  the  wave  length. 

For  the  subsequent  discussion  we  need  only  know  its  limiting 
behavior.  For  x  «  d,  R  1;  and  for  x  »  d,  R  — »  0.  As 
x  varies  from  x  «  d  to  x  »  d,  one  would  expect  a  monotonic 
transition  between  these  two  limiting  values  of  R. 

Consider  the  phase  changes  that  occur  when  x  is  increased.  For 
kx  «  1  and  R  »  1  we  have  approximately 

ps  -  upcQ/-ikx  . 

Thus,  the  phase  shift  between  u  and  p„  is  tt/2.  For  R  «=  0  we  have 

J 

P  =  upc  exp  ikx 

5 
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The  phase  shift  is  kx.  If  the  transition  from  R  «  1  to  R  «  0 
occurs  before  kx  =  tt/2,  there  must  be  a  minimum  in  the  phase 
shift  as  shown  in  Fig.  13.  This  phase  transition  is  indeed 
observed  (see  Fig.  8).  Note  that  the  point  where  R  0  is 
reached  depends  primarily  on  the  microphone  size. 

F .  A  Finite  Microphone  in  Front  of  an  Infinite  Transmitter 
in  a  Rarefied  Gas 

One  obtains  the  sound  propagation  by  multiplying  Eq.  (10.7) 
by  F(c)  and  integrating  over  c.  This  gives 


P ( s )  =  up  ^  exp[-(c/c  )2j 

>Ttt  J  o  c-  m 

m 


•  (1  +  H)  exp(iuye)  dc  _  (l0.8) 

1  -  R“  exp(2i^X/c) 

For  R  «  1  and  u*  «  c  we  can  approximate  this 

m 


.  r  •#  4 

P ( s )  --=  up  —  — 2- —  exp[-(c/c  )2]  dc 

Vtt^o  -ic~W  m 

m 


or 


P(s) 


•a  upc 

3  r  m 
2  -i^x 


UPC0 

-ikx 


(10.9) 


(10.10) 


Equation  (10. 10)  differs  from  the  corresponding  sound  wave  solu¬ 
tion  only  by  a  factor  9/5. 


-35- 


Report  No.  1169 


Bolt  Beranek  and  Newman  Inc 


0.  An  Attemrt  to  Explain  the  Measured  Results 

Calculated  curves  for  limiting  cases  are  presented  in  Pig.  14. 

In  order  to  normalize  the  data,  the  square  of  the  sound  pressure 

2  2  2 

at  each  separation  is  divided  by  u  p  cffl  ,  As  before,  s  * 

The  decay  from  s  ■  1  to  s  *  8  is  calculated  from 

|p(s)|2  *  !p(s  -  1)|2  exp  {-2  e0Axj 

-  | p ( s  -  1 ) | 2  exp  |-0.4vH57!>  (s  -  l)j  (10.12) 

by  using  the  value  a /$  *  0.2  (c.f.  Pig.  3).  The  plot  of  the 

result  of  this  calculation  is  labeled  Eq.  (10.12)  in  Fig.  14. 

For  the  decay  in  the  region  s  «  1  we  take  R  *  1  and  use 

Eq.  (10.10).  This  calculation  is  represented  by  the  portion 
of  the  curve  labeled  Eq.  (10,10).  Between  the  calculated  curves 
Eq.  (10.12)  and  Eq.  (10.10)  extrapolated  transitions  can  be  con¬ 
structed.  The  transition  curves  depend  on  the  transition  from 
R  ■  1  to  R  ■  0,  The  two  extreme  cases  of  transition  are  depicted 
in  Pig.  14,  that  of  a  fast  transition  from  R  *=  1  to  R  *  0  and 
that  of  8 low  transition  between  these  limits. 

Figure  15  shows  a  comparison  between  the  estimated  and  measured 
values.  The  measured  values  were  taken  from  Fig.  9.  The  results 

favor  a  fast  transition  from  R  -  1  at  s  «  1  to  R  ■  0  for  s  >  1. 

♦ 

Tne  pressures  at  x  *  0.05  cm  are  adjusted  to  be  equal.  Therefore, 
u/»  must  be  the  same  in  both  cases.  This  implies  that  in  the  10  kc 
case  the  amplitude  of  the  transmitter  is  five  times  as  high  as 
that  of  the  2  kc  case.  Thus,  one  would  expect  the  10  kc  curve  to 


-36- 


Report  No.  1169 


Bolt  Beranek  and  Newman  Inc 


be  approximately  14  db  higher  than  tne  2  kc  curve  for  s  >  1, 
i.e.,  when  the  Initial  decay  region  has  been  passed.  The 
difference  between  the  two  curves  at  x  =  3  cm  is  10  db.  We  add 
to  that  5  db.  This  represents  the  excess  attenuation  of  the  10  kc 
wave  over  the  2  kc  wave  in  reaching  the  point  x  «  3  (s  ■  5.1). 

Then  the  data  agree. 
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XI.  REFLECTION  OF  SOUND  FROM  SOLID  SURFACES 

The  results  of  experiments  for  both  the  amplitude  and  phase, 
expecially  at  small  separations,  suggest  the  possibility  of 
observing  reflection  of  sound  from  solid  surfaces  at  low  pres¬ 
sures.  Also  recall  that  the  experimental  results  favor  the 
conclusion  that  interaction  of  molecules  with  solid  surfaces 
is  predominantly  diffuse  elastic.  This  is  in  agreement  with 
the  finding  of  Stickney*^  who  used  a  molecular  beam  and  a  Rayleigh 
disc  to  measure  the  momentum  transfer  between  a  gas  and  metallic 
surface . 

We  set  up  an  experiment  to  examine  the  reflection  of  sound  di¬ 
rectly,  as  shown  in  schematic  form  in  Fig.  16.  With  this  arrange¬ 
ment  there  is  no  direct  path  between  the  transmitter  and  receiver. 
The  sound  reaches  the  microphone  primarily  by  reflection  from  the 
test  surface.  The  results  obtained  are  shown  in  Fig.  17.  The 
records  illustrate  clearly  that  reflection  of  sound  can  be  measured 
in  a  rarefied  gas.  Moreover,  the  results  in  Fig.  17  show  that 
standing  wave  patterns  of  a  sort  occur  also  at  low  pressures. 

The  maxima  and  minima  of  the  standing  waves  may  bo  used,  as  they 
are  used  at  higher  pressures,  to  determine  reflection  properties 
(e.g.,  relative  accommodation  coefficients)  of  surfaces  at  reduced 
pressures. 

If  one  assumes  that  the  distance  between  two  minima  is  directly 
proportional  to  c/u),  c  being  the  phase  velocity  of  propagation, 
then  the  results  in  Fig.  17  imply  that  the  phase  velocity  of 
sound  increases  with  a  decrease  in  pressure.  This  is  in  agree¬ 
ment  with  the  phase  variation  with  pressure  reported  in  Figs.  10 
and  11. 
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More  refined  techniques  would  be  required  to  make  the  measurements 
useful.  However  the  results  obtained  thus  far  do  give  encourage¬ 
ment  that  such  an  approach  may  be  successful. 
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XII.  VARIATION  OP  SOUND  PRESSURE  LEVEL  AS  A  FUNCTION  OF  PRESSURE 

The  reeulta  obtained  for  the  variation  of  sound  pressure  level 
as  a  function  of  separation  (e.g.,  see  Fig.  7)  indicate  that  the 
decay  at  small  separation  Is  Independent  of  the  pressure  (in  the 
range  of  multiple  reflections).  To  examine  this  phenomenon  in 
a  more  direct  way,  we  measured  the  variation  of  the  sound  pres¬ 
sure  level  as  a  function  of  pressure  at  a  frequency  of  6  kc  and 
a  constant  separation  of  0.35  cm  (s  *  0.32). 

The  result  of  this  experiment  is  shown  in  Fig.  18.  A  curve  based 
on  a  linear  relationship  is  also  plotted.  The  range  of  pressure 

_2i 

covered  by  this  experiment  extends  from  6  mm  Hg  down  to  1.4  x  10 

mm  Hg.  The  pressures  were  measured  by  an  ionization  gauge  from 
-4  .■* 

1.4  x  10  mm  Hg  to  about  5  x  10  J  mm  Hg  and  with  an  alphatron 

gauge  in  the  range  from  7  x  10“  3  mm  Hg  to  6  mm  Hg.  Note  that 

slight  discrepancies  between  theory  and  experiment  occur  at  the 

limit  of  the  gauge  range  (for  the  ionization  gauge  the  upper  limit 

is  at  5  x  10‘3  mm  Hg)  or  where  a  switch  of  scale  takes  place  (the 

alphatron  scale  has  to  be  switched  at  intervals  of  a  factor  of 

ten  in  pressure  starting  at  10  mm  Hg).  These  are  the  positions 

where  one  may  expect  inaccuracies  in  the  gauges  to  be  most  pro- 

-2 

nounced.  The  break  in  the  curve  at  about  10  mm  Hg  also  shows 
that  the  calibration  of  the  two  gauges  differs  by  a  small  amount. 
The  transition  region  from  the  classical  to  the  geometrical 
relaxation  regime  occurs  in  the  range  of  pressure  of  about  1.5  x 
10  mm  Hg.  Since  the  curve  does  not  exhibit  a  changing  slope, 
our  argument  that  no  change  would  be  expected  in  the  range  where 
s  is  much  smaller  than  unity  is  substantiated. 
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XIII.  EXPERIMENTS  POR  THE  MEASUREMENT  OF  SURFACE 
ACCOMMODATION  COEFFICIENTS 

A.  Outline  of  the  Theory 

A  configuration  useful  for  the  study  of  surface  reflections  Is 
given  in  Pig.  19.  is  the  surface  under  study.  S2  is  a  con¬ 
trol  surface  used  to  define  a  source  of  coherent  momentum  fluctua¬ 
tions  and  S^  is  a  control  surface  used  to  define  the  coherent 
part  of  the  distribution  function  of  those  molecules  reflected 
from  S^. 

As  in  the  previous  development,  we  ignore  the  source  terms  in 
Eq.  (1.1)  and  subject  the  homogeneous  equation  to  boundary  con¬ 
ditions  appropriate  to  the  configuration  under  study.  Then  the 
relationship  between  the  distribution  function  at  phase  point 
r  ,6  at  time  t  to  the  distribution  functions  at  the  phase  point 

r_L  at  t  in  a  region  having  no  collisions  may  be  described  by  a 
m 

propagator  T„fr£  t  Ir  it)  as  in  Eq.  (4.7).  n  denotes  the  trans- 
run Vn-n  n'-jiu-m  m 

mitting  control  surface  and  m  the  receiving  control  surface.  However, 
when  studying  surface  interactions,  one  must  give  explicit  considers- 

t 

tion  to  the  transformation  caused  by  the  surface  itself.  By  denoting 

the  distribution  function  of  those  molecules  striking  surface  by 

P  q  and  those  leaving  by  F  1  Q  ,  one  can  indicate  the  surface  trans- 
-bl  ^1 

formation  symbolically  by 

P+Si  *  *  (13,1) 


It  is  to  the  determination  of  MQ 

bl 


thi3  section. 


that  we  direct  our  efforts  in 
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The  control  surface  Sg  is  essentially  a  source  of  "test"  molecules, 
l.e.  molecules  having  a  coherent  momentum  fluctuation  characterized 
by  the  fluctuating  velocity  a  .  The  control  surface  is  estab¬ 
lished  to  assess  the  distribution  function  and/or  velocity  moments 
of  those  molecules  arriving  from  which  have  retained  this 
coherent  momentum.  Thus  it  is  necessary  to  express  P  „  ,  the  dis- 
tribution  of  molecules  arriving  at  in  terms  of  Jvia  the 

propagators  T  and  the  surface  transformation  M<,  .  This  relation¬ 

al 


ship  is  expressed  as  follows: 


V  ‘/•••/{dT}Tl*1T21F+S2 


•  (£2,  ^2,  0,  t2)  .  (13.2) 

M  represents  the  elemental  variables  of  all  initial  and  inter¬ 
mediate  phase  points  on  the  appropriate  surfaces.  As  in  Eq.  (3.1) 

P  „  can  be  separated  into  a  dominant  incoherent  part  F°  and  a 
~S3 

coherent  perturbation  h.  Thus, 


f_s3  ^  $3* 


(13.3) 


We  first  determine  whether  measurements  of  h  or  of  moments  over  h 
can  determine  parameters  which  can  be  used  to  characterize  M,  the 
surface  interaction  operator. 


It  is  convenient  to  express  the  operator  M  as  a  sum  of  three  terms 


M 


»s  ♦  Mde  +  Mdl 


(13.4) 
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Each  term  is  related  to  one  of  the  three  processes  of  surface 
interaction  as  discussed  in  Section  IV. C.  Expressing  M  in  this 
manner  amounts  to  the  assumption  that  the  three  processes, 
specular,  diffuse  elastic,  and  diffuse  inelastic  reflections, 
are  linearly  independent.  Such  an  assumption  may  be  valid  pro¬ 
vided  one  seeks  to  describe  only  the  gross  properties  of  the 
molecular-surface  interactions.  Mg,  Mde  and  Mdi  are  the  operators 
describing  the  processes  of  specular,  diffuse  elastic  and  diffuse 
inelastic  interactions  respectively. 


Equation  (13.*0  can  be  substituted  in  Eq.  (13.2)  to  yield  a 
decomposition  of  h,  Eq.  (13.3),  into  three  parts: 


h 


hs  +  hde  +  h 


‘di 


(13.5) 


The  receiver  is  in  general  a  momentum  sensitive  device  and  i3 

capable  therefore  only  of  measuring  a  linear  combination  of  moments 

of  the  distribution  function.  The  suitability  of  an  experiment  to 

determine  surface  interaction  parameters  is  dependent  on  the  manner 

in  which  the  properties  of  influencr  the  velocity  moments  at 

51 

S3.  The  result  is  a  function  of  2111(1  ®13  as  evidenced  by  the 
expression  Eq.  (13.2).  Consider  as  a  representative  moment,  the 
momentum  flux  calculated  from  Eq.  (13.5).  Since  each  of  the  con¬ 
tributing  components  of  h  is  a  different  function  of  angle,  the 
total  momentum  can  be  written 

pS3  =  RsQs^21,U13^  +  Rde^de^021,913^  ’  (13.6) 

Rg  and  Rde  are  the  measures  of  the  specular  and  diffuse  elastic 
reflection  and  are  related  to  the  momentum  accommodation  coefficient 
of  the  two  processes.  ftg  and  are  in  general  complex  functions 


-43- 


Report  No.  1169 


Bolt  Beranek  and  Newman  Inc 


of  the  sizes  of  surfaces  S2  and  S^.  The  tern  Involving  the  dif¬ 
fuse  inelastic  reflection  is  missing  in  Eq.  (13.6)  since  this 
process  destroys  the  coherent  fluctuating  momentum  of  the  imping¬ 
ing  molecules. 

Equation  (13.6)  is  of  exceptional  complexity.  It  is  therefore 
necessary  to  introduce  constraints  related  to  specific  experiments 
so  designed  that  their  interpretation  can  be  related  to  a  simpli¬ 
fied  form  of  Eq.  (13.6).  Such  simplification  can  be  obtained  if 
the  surfaces  S2,  and  can  be  made  small  enough  30  that  the 
path  of  propagation  from  S2  to  can  be  specified  to  a  good 
approximation  by  ®21  and  only.  A  schematic  representation  of 
the  experimental  system  is  shown  in  Pig.  20a.  For  such  a  system, 
Eq.  (13.6)  can  be  approximated 

Ps3  *  ^3^8^21' ^13)  +  Rdecos ^13^  •  (!3.7) 

Here  ^g(02l,0i3^  13  a  func^lon  such  that  it  is  equal  to  unit 
when  92i  m  ^13  3111(1  approaches  zero  rapidly  as  |®21  - 
increases.  An  approximate  expression  for  thi3  function  can  be 
formally  written 


Q  ) 
2V  13' 


sm(e2i  - 


J21  ■  ®13 


2n 


(13.8) 


where  n  is  a  positive  integer.  ®21  and  in  Eq.  (13.8)  are 
expressed  in  radians.  In  Eq.  (13.7)  we  assume  that  diffraction 
effects  from  the  edges  of  surface  are  negligible. 
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By  keeping  both  the  distances  and  the  source  output  constant, 
measurements  of  pc^  as  a  function  of  the  angle  should  lead 
to  determination  of  such  ratios  as 


pS3(tf13  “  *12)  Rs  c°3 ( ©12 ) 

’  T^  +  cosT^7 


(13.9) 


from  which  the  ratio  Rg/Rde  can  be  ascertained.  If  a  baffled 
opening  is  used  to  replace  the  reflector  surface,  an  arrangement 
that  is  illustrated  in  Fig.  20b,  the  measured  sound  pressure  at 
is  proportional  to 


(e; 


13 


V 


Rs  +  Rde  +  Rdl 


(13.10) 


When  the  source  strength  and  distances  involved  are  maintained 
equal  in  both  experiments,  an  additional  ratio  can  be  obtained, 
namely 


PS. 


'1 


(°13  ** 


12 


Uiprrj  -  "de 


cos(013) 


(13.11) 


Hence,  one  should  obtain  the  relative  mangitudes  of  R 


,  Rde  and 


di* 


The  baffled  hole  experiment  can  serve  further  as  a  check 


that  diffraction  effects  are  indeed  negligible.  This  is  achieved 
by  measuring  ensurln8  that  p^  (0^3  /  °12^^ 

P33^13  =  J12^’ 
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B.  Outline  of  the  Experiment 

The  schematic  arrangement  given  In  Fig.  20  shows  the  experimental 
setups  Sp  represents  a  1.25-inch  transmitter,  represents  a 
1-inch  circular  aluminum  disk,  Sj  represents  a  1-inch  hole  in  the 
center  of  a  5-inch  aluminum  baffle  and  represents  a  1-inch 
microphone.  The  distances  and  ^3  are  each  3-lnches.  Is 
attached  to  a  3-inch  long  arm  (1^)  that  can  be  actuated  externally 
when  the  system  is  In  the  bell  jar.  The  arm  can  be  actuated  to 
move  the  microphone  In  a  circular  path  In  a  plane  containing  the 
centers  of  S2*  S1  (or  S£)  and  (or  Sp  with  S1  as  the  center 
point.  The  geometrical  relationship  between  the  various  com¬ 
ponents  are  thus  not  fine  enough  to  make  Eqs.  (13.7)  and  (13.10) 
applicable  In  an  absolute  sense,  nevertheless  we  hoped  that  a 
rough  approximation  could  be  achieved.  However,  the  signal- to- 
nolse  ratio  was  too  low  and  no  dependable  results  are  reported. 

To  Improve  the  signal- to-noise  a  more  powerful  transmitter  was 
designed  and  constructed.  With  the  more  powerful  transmitter, 
the  signal -to-nolse  ratio  is  still  too  low.  It  appears  that  at 
least  10  to  20  db  higher  signals  are  required  In  order  to  obtain 
meaningful  results.  This  Is  not  possible  with  our  existing  tools. 
We  believe  there  is  yet  the  possibility  of  carrying  out  such  an 
experiment  successfully.  Two  suggestions  follow: 

1.  Replace  the  source  at  Sg  by  a  reasonably  collimated  and 
powerful  molecular  beam  and  place  a  transmitter  at  S^. 
Equation  (13*6)  should  then  be  modified  to  accommodate 
the  new  arrangement.  Such  a  modification  can  be  derived 
using  the  methods  of  Section  IV. B. 
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2.  Replace  the  source  at  Sg  by  a  collimated  and  powerful 
molecular  beam  that  Is  chopped  at  a  fixed  frequency. 

By  the  time  the  fluctuating  disturbance  reaches  the 
surface  the  large  amplitudes  would  have  diminished 
sufficiently  to  be  considered  small  and  the  molecules 
will  have  predominantly  a  fluctuating  component  of 
momentum  at  the  frequency  of  the  chopper.  Both  the 
colllmatlon  and  the  higher  amplitudes  that  can  be 
obtained  at  by  this  means  should  contribute  more 
than  20  db  to  the  signal  to  noise  ratio  making  the 
experiment  feasible. 
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XIV.  A  KINETIC  THEORY  MODEL  OP  SOUND  PROPAGATION 

There  is  a  close  quantitative  correspondence  between  the  results 
of  Section  VI  which  we  obtained  for  the  geometrical  relaxation 
regime  and  those  which  are  obtained  in  the  classical  and  frequency 
relaxation  regimes. The  point  of  departure  i3  that  in  the  lat¬ 
ter  regimes  one  must  identify  the  variable  s  with  the  mean  free 
path,  s  ■  caA/cm,  whereas  in  the  geometrical  relaxation  regime  it 
is  defined  in  terms  of  the  normal  distance  from  the  transmitter  to 
the  point  of  observation,  Eq.  (5.2).  However,  in  both  cases  these 
lengths  represent  the  free  flight  path  between  collisions.  Since 
the  dispersive  properties  of  sound  in  the  geometrical  relaxation 
regime  are  related  to  the  free  flight  path,  we  infer  from  this 
correspondence  that  dispersion  in  the  other  regimes  can  be  related 
to  the  free  flight.  We  propose  the  following  model.  In  the 
classical  and  frequency  relaxation  regimes  a  coherent  ensemble 
of  molecules  travels  an  average  distance  of  one  mean  free  path 
without  collision.  In  a  collision  process  a  new  ensemble  of 
molecules  is  statistically  established  to  which  the  coherence  is 
communicated.  Thl3  ensemble  of  molecules  now  becomes  the  carrier 
of  this  portion  of  the  3ound  field.  Between  collisions,  the 
motion  is  essentially  in  the  geometrical  relaxation  regime  and 
is  subjected  to  the  loss  of  coherence  typical  of  that  regime. 

In  a  subsequent  collision  process  the  residual  coherence  is  once 
more  established  in  an  ensemble  of  molecules  which  has  a  basically 
similar  molecular  distribution  to  the  preceding  ensembles;  and  so 
on.  Thus,  with  the  proper  interpretation  as  to  what  constitutes 
sound  propagation  in  the  classical  and  the  frequency  relaxation 
regimes,  a  3imple  model  can  be  constructed  to  account  for  the 
dispersive  properties  of  sound  in  these  regimes. 
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A  simple  way  of  constructing  the  model  Is  to  envision  collision 
control  surfaces  as  depicted  in  Fig.  21.  We  characterize  the 
collision  process  and  the  subsequent  emergence  of  a  fluctuating 
component  of  velocity  by  a  function  of  velocity  0(v  ).  At  the 
zeroth  control  surface  the  fluctuating  component  in  the  density 
pertaining  to  those  molecules  that  travel  in  the  positive  x- 
direction  is  given  by 

P O*/  0<vx>  F<vx>  dvx  *  f14*1) 

o 

where  F(vx)  is  the  distribution  function  of  molecular  velocities 
in  the  positive  x-direction.  In  travelling  one  mean  free  path 
from  the  zeroth  to  the  first  control  surface,  the  corresponding 
density  fluctuation  at  the  first  control  surface  is  given  by 

Pi  -  f  0(vx)  F(vx)exp(iAo>/vx)  dvx  .  (14.2) 

o 

Proceeding  in  this  manner,  it  can  be  readily  shown  that  the  den¬ 
sity  fluctuation  at  the  Nth  surface  is  given  by 

PN  -  PotPi/Pol”  •  (I1*. 3) 

The  fluctuating  component  in  the  distribution  function  leaving 
the  Nth  surface  is  given  by 

hN  "  tPi/P0)N  0(vx>  P(vx>  ' 
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Denote  the  distance  from  the  Initial  control  surface  where  the 
sound  Is  generated  to  a  control  surface  where  the  sound  Is 
finally  measured  by  x.  Then 

x  -  NA  +  »x  ,  (14.5) 

where  &x/A  <  1.  Since  the  sound  measuring  device  Is  a  momentum 
flux  (pressure)  assessor,  the  relationship  between  the  sound 
pressure  at  x  -  0,  pQ,  and  x  ■  N  +  Ox,  p(x),  Is  given  by 

p(x)  -  P0lPi/P</  /  *<vx>  vl  p<vx> 

r  0  o 

•  exp(l&x*J>/vx)  dvx  .  (14.6) 

If  N  Is  large,  Eq.  (14.6)  can  be  readily  approximated  by 

p(x)«  P0tP1/Pol (jc/A)  .  (14.7) 

Thus  In  the  frequency  relaxation  regime  the  pressure  and  the 
density  fluctuations  obey  essentially  the  same  dispersive  laws. 
This  Is  of  course  In  sharp  contrast  to  the  geometrical  relaxa¬ 
tion  regime  where  the  density  fluctuations  and  the  pressure 
fluctuations  do  not  necessarily  obey  the  same  dispersive  laws. 

If,  in  an  experiment,  the  planes  perpendicular  to  the  direction 
of  propagation  where  observations  of  the  sound  pressure  are  made 
are  separated  such  that  Is  large, 

N21  ■  —  A  *X  »  1  »  '  (n*.8) 
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where  x1  la  the  location  defining  the  first  plane  and  x2  the 
second  plane ,  then  the  attenuation  and  phase  parameters  are  given 
by 


°  -  IP  -  -  j  ^n(p1/p0) 


(1H.9) 


Alternately  If  N21  <  1,  the  attenuation  and  phase  parameters 
are  similar  to  those  of  the  geometrical  relaxation  regime,  namely. 


a  -  ip 


/  ^x)vxp(vx>«p("Ax)  dvx 
m/o  0(vx)v^P(vx)exp(lAxa>/vx)  dvx 


(14.10) 


where  Ax  «  xg  -  x1  .  Observe  that  the  attenuation  and  phase 
parameters  In  Eq.  ( 14.10)  are  independent  of  the  ambient  pres¬ 
sure. 


For  the  transition  region  from  N21  »  1  to  N21  **  1  the  disper¬ 
sion  properties  can  be  derived  from  Eq.  (14.6);  however,  the 
expressions  in  this  case  are  more  complex. 

To  calculate  further  one  must  explicitly  know  the  function 
0(v  )  .  Thi3  function  embodies  the  collision  processes  of 
molecular-molecular  interactions.  In  Section  VI  we  showed  that 
the  dispersive  properties  are  not  very  sensitive  to  the  details 
of  the  dependence  of  the  interactions  on  the  velocity  vx,  and 
consequently  one  would  expect  Eq.  (14.9)  to  yield  results  that 
are  not  very  sensitive  to  the  explicit  dependence  of  0(vx)  on 
v  .  This  proves  to  be  true.  For  any  finite  polynomial  function 

4^ 
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4(vx),  Eq.  (14.9)  yields  attenuation  and  phase  parameters  that 
show  similar  behavior  to  that  of  the  corresponding  geometrical 
relaxation  parameters  depicted  in  Fig.  3#  with  A  in  place  of  x. 

The  experimental  results  support  this  conclusion,  as  is  apparent 
from  the  data  of  Figs.  3*  10  and  11. 

Although  there  are  theories  which  lead  to  dispersion  laws  for 
normal  gases  and  plasmas.  It  remains  a  point  of  continued  dis¬ 
cussion  to  obtain  a  physical  description  of  the  dispersive 
properties  uncovered  by  these  formal  theories.  The  procedure 
that  we  have  suggested  above  has  the  advantage  that  the  physical 
Interpretation  Is  present  In  the  model.  We  emphasize  the  essential 
distinction  between  a  model  of  this  type  which  Is  motivated  by  its 
physical  content  and  other  models  that  have  been  used  in  this  type 
of  problem  which  are  motivated  more  by  their  solvability.  The 
conceptual  basis  Is  simply  that  the  dispersive  properties  of 
sound  In  gases  is  directly  determined  from  an  assessment  of  the 
average  amount  of  ’’free  flight”  that  the  molecules  experience. 
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FIG.  1  THE  3  REGIMES  OF  SOUND  PROPAGATION 
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FIG.  3  COMPUTED  ATTENUATION  AND  PHASE  FOR  SEVERAL  TYPES 
OF  SURFACE  INTERACTIONS  (o  EXPERIMENTAL  VALUES. 


FIG.4  SCHEMATIC  SKETCH  OF  THE  PUMPING  SYSTEM 
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FIG. 7  TYPICAL  RECORDS  OF  SOUND  PRESSURE  LEVEL 
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FIG. 8  TYPICAL  RECORDS  OF  PHASE  VS.  SEPARATION 
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FIG. 9  INITIAL  DECAY  CURVES  FOR  DIFFERENT 
FREQUENCIES  (10  kc  AND  2  kc)  AT 
PRESSURES  OF  0.  46  m 
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ATTENUATION  AND  PHASE  WITH  PRESSUR 
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FIG.  12  TRANSMITTER  -  MICROPHONE  SYSTEM 
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FIG.  14  RELATIVE  SOUND  PRESSURE  LEVEL  AS 
A  FUNCTION  OF  THE  NORMALIZED 

SEPARATION  PARAMETER  $ 
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FIG. 15  RELATIVE  SOUND  PRESSURE  LEVEL  AS 
A  FUNCTION  OF  THE  NORMALIZED 
PARAMETER  S 
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FIG.  16  TRANSMITTER  -  REFLECTING  SURFACE - 
MICROPHONE  SYSTEM 


RELATIVE  SOUND  PRESSURE  LEVEL 


REPORT  NO.  1169 


BOLT  BERANEK  ft  NEWMAN  INC 


COATED 

MYLAR 


COATED  MYLAR 
AND  GAUZE 


-v/  w 


A: 


COATED  MYLAR 
AND  GAUZE 


COATED 

MYLAR 


10  DB 


6.5 


7.5 


i 

8.5cm 


FIG.  17  REFLECTION  OF  SOUND  FROM  SURFACES  AS 
A  FUNCTION  OF  SEPARATION 
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FIG. 18  SOUND  PRESSURE  LEVEL  VS.  PRESSURE 

(FREQUENCY  -  6  kc,  SEPARATION  *  0.35  cm, 

MICROPHONE  DIAMATER  *  1  IN.,  TRANSMITTER 
DIAMETER  *  5.0  IN.) 
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FIG.  19  REFLECTION  OF  SOUND.  S,  A  TRANSMITTER 
SURFACE,  S,  A  REFLECTING  SURFACE  AND 
S}  A  RECEIVER  SURFACE 


FIG.  20  SCHEMATIC  REPRESENTATION  OF  AN 

APPARATUS  FOR  STUDIES  OF  MOLECULAR- 
SURFACE  INTERACTIONS 
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